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Abstract 

Armchair graphene nanoribbon (AGNR) for n = 3m and n = 3m + 1 family carrier statistic under uniaxial strain is 
studied by means of an analytical model based on tight binding approximation. The uniaxial strain of AGNR carrier 
statistic models includes the density of state, carrier concentration, and carrier velocity. From the simulation, it is found 
that AGNR carrier concentration has not been influenced by the uniaxial strain at low normalized Fermi energy for 
n = 3m and n = 3m + 1 . In addition, the carrier velocity of AGNR is mostly affected by strain at high concentration of 
n & 3.0 x 1 0 7 and 1 .0 x 1 0 7 m _1 for n = 3m and n = 3m + 1 , respectively. The result obtained gives physical insight 
into the understanding of uniaxial strain in AGNR. 
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Background 

Graphene has attracted numerous research attention 
since it was isolated in 2004 by Novoselov et al. [1]. Due 
to its unique hexagonal symmetry, graphene posses many 
remarkable electrical and physical properties desirable in 
electronic devices. It is the nature of graphene that it does 
not have a bandgap, which has limited its usage. There- 
fore, efforts to open up a bandgap has been done by several 
methods [2-4]. The most widely implemented method 
is patterning the graphene into a narrow ribbon called 
graphene nanoribbon (GNR) [4]. Recently, strain engi- 
neering have started to emerge in graphene electronics 
[5] . It is found that strain applied to graphene can modify 
its band structure, thus, altering its electronic properties 
[6-8]. In fact, uniaxial strain also helps in improving the 
graphene devices electrical performance [9] . Similar char- 
acteristics have been observed when strain is applied to 
conventional materials like silicon (Si), germanium (Ge), 
and silicon germanium (SiGe) [10]. Strain in graphene can 
be characterized by two major varieties, namely uniax- 
ial and shear. This strain behaves differently on graphene 
depending on the edge shape, namely zigzag or arm- 
chair [8] . The presence of the strain effect in graphene is 
by the G peak that splits and shifts in the Raman spec- 
trum [11,12]. It is worth noting that strain in graphene 
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may unintentionally be induced during the fabrication of 
graphene devices. 

Computational modeling and simulation study pertain- 
ing to strain graphene and GNR for both the physical and 
electrical properties have been done using few approaches 
such as the tight binding model and the ab initio calcu- 
lation [6,13]. An analytical modeling approach has also 
been implemented to investigate the strain effect on GNR 
around the low-energy limit region [14,15]. However, 
most of the previous works have only focused on the elec- 
tronic band structure, particularly the bandgap. As the 
carrier transport in GNR has a strong relation with this 
electronic band structure and bandgap, it is mandatory to 
investigate the strain effect on the carrier transport such 
as carrier density and velocity. Therefore, in this paper, 
an analytical model representing uniaxial strain GNR car- 
rier statistic is derived based on the energy band structure 
established by Mei et al. [15]. The strain effect in our 
model is limited to low strain, and only the first subband 
of the AGNR n = 3m and n = 3m + 1 families is consid- 
ered. In the following section, the analytical modeling of 
the uniaxial strain AGNR model is presented. 

Methods 

Uniaxial strain AGNR model 

The energy dispersion relation of GNR under tight bind- 
ing (TB) approximation incorporating uniaxial strain is 
represented by Equation 1 taken from reference [15]. The 
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Figure 1 Energy band structure of uniaxial strain AGNR (a) n = 3m and (b) n = 3m + 1 for the model in Equation 1 . 



TB approximation is found to be sufficient in the investi- 
gation for small uniaxial strain strength. This is because 
the state near the Fermi level is still determined by the 2p z 
orbitals that form the n bands when the lattice constant 
changes [6]: 
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where h = -^ £i t 2 = t 3 = p^jj, to = -274 eV is the 
unstrained hopping parameter, a — 0.142 nm is the lat- 
tice constant and t\ and t 2 are the deformed lattice vector 
hopping parameter of the strained AGNR. s is the uniaxial 
strain [15]. 

Using the first-order trigonometric function, Equation 1 
can further be simplified to the following equation: 
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To model the bandgap, at k x = 0, Equation 2 is reduced 
to [15] 
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Thus, the bandgap is obtained as the following equation 
[15]: 



E g = 2E (k x = 0) = 2 
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The energy dispersion relation from Equation 2 can 
further be simplified to 



E(n,k) = ± 
where 
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Equation 5 will be the basis in the modeling of strain 
GNR carrier statistic. GNR density of state (DOS) is fur- 
ther derived. The DOS that determines the number of 
carriers that can be occupied in a state of the system [16] 



(a). 



<s 2 

C 

PQ 1 



unstrained 

8=3% 

Si(1.12eV) 
Ge (0.67eV) 



10 



GNR Width (nm) 



(b) 



— unstrained 
—8=2% 

8=4% 

— 8=6% 




2 4 6 

GNR Width (nm) 



Figure 2 Bandgap of AGNR in respond to the width for (a) n = 3m and (b) n = 3m + 1. 



Johari and Ismail Nanoscale Research Letters 201 3, 8:479 
http://www.nanoscalereslett.eom/content/8/1/479 



Page 3 of 6 



(a) 



x 10 



> 4 

O 9 
Q 2 



0 



-0.5 



. unstrained 

■ 8=1% 

■ 8=2% 

■ 8=3% 



0 



0.5 



(b) 



x 10 

— unstrained— 8=2% 8=6% 
- 8=4% 



> 2 

w 

go 

° 1 
Q 1 



0 
-0.5 



l_U 



0 



0.5 



Energy (eV) Energy (eV) 

Figure 3 DOS varying the uniaxial strain strength in AGNR (a) n = 3m and (b) n = 3m + 1. 



is yielded as in Equation 7: 



DOS(£) 



(7) 



E 2 



In the modeling of the strain GNR carrier concen- 
tration, energy dispersion relation is approximated with 
the parabolic relation, \fl + a 2 & +\a. By substitut- 



ing the normalized Fermi energy as rj 



E F -E g /2 



the 



strain AGNR carrier concentration model is derived and 
represented by 



(8) 



To further evaluate the intrinsic carrier velocity in 
response to the uniaxial strain, the following definition is 
referenced [17]: 



v 7 /DOS(£) xf(E) 
DOS(£) xf(E) 
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The Fermi velocity, vj, is modeled as in reference [18]. 
Thus, Vf is obtained as the following equation: 



E 2 



(10) 



Hence, using the intrinsic velocity model defined in 
Equation 9, the strain AGNR intrinsic carrier velocity 
yields the following equation: 



2n^BVhT goW 



(ID 



The analytical model presented in this section is plotted 
and discussed in the following section. 

Results and discussion 

The energy band structure in respond to the Bloch wave 
vector, k x , modeled as in Equation 1 which was established 
by Mei et al. [15], is plotted in Figure 1 for n = 3m and 
n = 3m + 1 family, respectively. For each simulation, only 
low strain is tested since it is possible to obtain experimen- 
tally [12]. It can be observed from both figures that there is 
a distinct behavior between the two families. For n = 3m, 
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Figure 5 Uniaxial strained AGNR carrier velocity in response to carrier concentration for (a) n = 3m and (b) n = 3m + 1. 



the separation between the conduction and valence bands, 
which is also known as bandgap, increases with the incre- 
ment of uniaxial strain. On the contrary, the n = 3m + 1 
family exhibits decrements in the separation of the two 
bands. It is worth noting that the n = 3m + 1 family also 
shows a phase metal-semiconductor transition where at 
7% of strain strength, the separation of the conduction and 
valence bands almost crosses at the Dirac point. This is 
not observed in the n = 3m family [15]. 

The hopping integral to between the n orbitals of AGNR 
is altered upon strain. This causes the up and down shift, 
the a* band, to the Fermi level, Ep [19]. These two phe- 
nomena are responsible for the bandgap variation. It has 
been demonstrated that GNR bandgap effect with strain 
is in a zigzag pattern [14]. This observation can be under- 
stood by the shifting of the Dirac point perpendicular to 
the allowed k lines in the graphene band structure and 
makes some bands closer to the Fermi level [7,8]. Hence, 
the energy gap reaches its maximum when the Dirac point 
lies in between the two neighboring k lines. The allowed 
k lines of the two families of the AGNR have different 
crossing situations at the K point [8]. This may explain 
the different behaviors observed between n = 3m and 
n — 3m + 1 family. 



To further evaluate, the GNR bandgap versus the GNR 
width is plotted in Figure 2. Within the uniaxial strain 
strength investigated, the bandgap of the n = 3m family 
is inversely proportional to the GNR width. The narrow 
bandgap at the wider GNR width is due to the weaker 
confinement [20]. The conventional material of Si and Ge 
bandgaps are also plotted in Figure 2 for comparison. In 
order to achieve the amount of bandgap similar to that of 
Si (1.12 eV) or Ge (0.67 eV), the uniaxial strain is projected 
to approximately 3% for the n = 3m family. A similar 
observation can be seen for n = 3m + 1 with 2% uniax- 
ial strain. However, a higher strain resulted in a different 
kind of observation. For example at 4% uniaxial strain, the 
phase transition from metallic to semiconductor occurs at 
a GNR width of approximately 3m, The phase transition is 
not observed in AGNR n = 3m [15]. When higher strain is 
applied, the phase transition occurs at a lower width. The 
difference in GNR width for the phase transition to occur 
depends on the subband spacing effect with GNR width 
[21]. The constitution of the phase transition suggests that 
the GNR bandgap can be tuned continuously between the 
metal and semiconductor by applying strain. 

Based on the energy band structure, the analytical 
model representing the DOS of strained AGNR is derived 




Fermi Velocity, ^ (ms ) x 10 
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Figure 6 Fermi velocity effect to the energy band structure of uniaxial strain AGNR for (a) n = 3m and (b) n = 3m + 1. 
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as in Equation 7. It is necessary to understand the DOS 
of strain AGNR as it will give insight on the amount of 
carriers that can be occupied in a state. The analytical 
model for strained AGNR is shown in Figure 3 for the first 
subband for the two AGNR families. It appears that the 
patterns of DOS are essentially the same for both AGNR 
families. It can be observed from Figure 3a,b that the 
Van Hove singularities are present at the band edge. For 
AGNR with n — 3m, the increment of strain increases 
the DOS remarkably. However, when s = 3%, despite the 
wide bandgap, the DOS substantially decreases. This is 
the reason for changing the band index, p, which corre- 
sponds to the bandgap [15]. In the case of n = 3m + 1, 
the DOS exhibits the opposite. In fact, when the strain 
strength made the band approach the transition phase, the 
DOS reduces significantly; at the same time, the bandgap 
approaches zero. 

To assess the effect of strain on AGNR carrier concen- 
tration, the computed model as in Equation 8 as a function 
of rj is shown in Figure 4. Apparently, the amount of car- 
riers increases when the AGNR n — 3m is added with 
uniaxial strain. Conversely, AGNR n — 3m + 1 shows 
a reduction in carrier concentration upon strain. Most 
notably, for AGNR n = 3m, the carrier concentration 
converges at low rj within the investigated strain level. 
Meanwhile, the carrier concentration exhibits consider- 
able effect upon the strain when the Fermi level lies at 
3I<bT away from the conduction or valence band edge. The 
same observation was achieve in AGNR n = 3m + 1. 

To assess the carrier velocity effect with carrier con- 
centration upon the strained AGNR, the analytical model 
in Equation 10 is plotted in Figure 5. It can be seen 
from Figure 5a,b that the GNR carrier velocity decreases 
and increases with the applied uniaxial strain for AGNR 
n = 3m and AGNR n = 3m + 1 families, respectively. 
Inspection of these figures also showed that the uniaxial 
strain mostly affected the carriers at high concentration. 
This is evident by the curves that tend to converge until 
n ^ 3 x 10 7 m -1 and has an almost constant veloc- 
ity at 1.8 x 10 5 ms -1 . When the concentration is high 
enough, the uniaxial strain starts to give a considerable 
effect to the velocity. This is supported by the previous 
observation in Figure 4 where the effect of the strain is 
infinitesimal at low 77. In fact, the applied strain also affects 
the degeneracy approach. The strained AGNR n = 3m 
approach degenerated later compared to the unstrained 
AGNR. A similar behavior was also observed in the AGNR 
n = 3m + 1 family except that strained AGNR approaches 
degeneracy faster compared to their unstrained counter- 
parts. This indicates that uniaxial strain is beneficial at a 
high concentration regime. Nonetheless, this is not unrea- 
sonable for low-dimensional nanostructures like GNR 
since it is mostly in the degenerated realm particularly for 
narrow width. 



The energy in response to the Fermi velocity of strained 
AGNR is shown in Figure 6. It can be observed that the 
effect of the strain on the Fermi velocity for both AGNR 
families is dramatic. Both AGNR n — 3m and n — 3m + 1 
have appreciable reduction in the Fermi velocity when the 
uniaxial strain increases as can be seen in Figure 6a,b. This 
reduction is attributed to the decrements in the n orbital 
overlap [22] in the AGNR band structure. As a conse- 
quence, the mobility is predicted to be degraded [23] as a 
result of the strong effect in the interaction of the strained 
carbon atoms [18,23]. 

Conclusions 

In this paper, the uniaxial strain AGNR for n = 3m and 
n = 3m + 1 family carrier statistic is analytically modeled, 
and their behaviors are studied. It is found that uniaxial 
strain gives a substantial effect to AGNR carrier statistic 
within the two AGNR families. The AGNR carrier con- 
centration has not been influenced by the uniaxial strain 
at low normalized Fermi energy. It is also shown that the 
uniaxial strain mostly affects carrier velocity at a high con- 
centration of n & 3.0 x 10 7 m -1 and n & 1.0 x 10 7 m -1 
for n — 3m and n = 3m + 1, respectively. In addition, 
the Fermi velocity of the AGNR n = 3m and n = 3m + 
1 exhibits decrements upon the strain. Results obtained 
give physical insight on the understanding of the uniaxial 
strain effect on AGNR. The developed model in this paper 
representing uniaxial strain AGNR carrier statistic can 
be used to further derive the current-voltage characteris- 
tic. This computational work will stimulate experimental 
efforts to confirm the finding. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

ZJ carried out the analytical modelling and simulation studies. Rl participated 
in drafting and improving the manuscript. Both authors read and approved 
the final manuscript. 

Acknowledgements 

The authors would like to acknowledge the financial support from the 
Research University grant of the Ministry of Higher Education (MOHE), Malaysia 
under project number R.J1 30000.7823.4F146 . Also, thanks to the Research 
Management Centre (RMC) of Universiti Teknologi Malaysia (UTM) for 
providing excellent research environment in which to complete this work. 

Received: 1 8 July 201 3 Accepted: 1 1 October 201 3 
Published: 14 November 2013 

References 

1 . Novoselov KS, Geim AK, Morozov SV, Jiang D, Zhang Y, Dubonos SV, 
Grigorieva IV, Firsov AA: Electric field effect in atomically thin carbon 
films. Science 2004, 306(5696):666-669. 

2. Castro EV, Novoselov KS, Morozov SV, Peres NMR, dos Santos JMBL, 
Nilsson J, Guinea F, Geim AK, Neto AHC: Biased bilayer graphene: 
semiconductor with a gap tunable by the electric field effect. Phys 
Rev Lett 2007, 99:216802. 

3. Nourbakhsh A, Cantoro M, Vosch T, Pourtois G, Clemente F, van der Veen 
MH, Hofkens J, Heyns MM, Gendt SD, Sels BF: Bandgap opening in 



Johari and Ismail Nanoscale Research Letters 201 3, 8:479 
http://www.nanoscalereslett.eom/content/8/1/479 



Page 6 of 6 



oxygen plasma-treated graphene. Nanotechnology 2010, 
21(43):435203. 

4. Li X, Wang X, Zhang L, Lee S, Dai H: Chemically derived, ultrasmooth 
graphene nanoribbon semiconductors. Science 2008, 31 9:1 229-1 232. 

5. Pereira VM, Neto AHC: Strain engineering of graphene's electronic 
structure. Phys Rev Lett 2009, 103(4):046 801 +. 

6. Gui G, Li J, Zhong J: Band structure engineering of graphene by 
strain: first-principles calculations. Phys RevB 2008, 78(7) 075435. 

7. Rosenkranz N, Mohr M, Thomsen C: Uniaxial strain in graphene and 
armchair graphene nanoribbons: an ab initio study. Annalen der 
Physik 201 1, 523(1 -2): 1 37- 1 44. 

8. Li Y, Jiang X, Liu Z, Liu Z: Strain effects in graphene and graphene 
nanoribbons: the underlying mechanism. Nano Res 201 0, 
3(8):545-556. 

9. Alam K: Uniaxial strain effects on the performance of a ballistic top 
gate graphene nanoribbon on insulator transistor. Nanotechnol 
IEEE Trans 2009, 8(4):528-534. 

1 0. Lee ML, Fitzgerald EA, Bulsara MT, Currie MT, Lochtefeld A: Strained Si, 
SiGe, and Ge channels for high-mobility metal-oxide-semiconductor 
field-effect transistors. JAppI Phys 2005, 97(1 ):01 1101. 

1 1 . Mohiuddin TMG, Lombardo A, Nair RR, Bonetti A, Savini G, Jalil R, Bonini N, 
Basko DM, Galiotis C, Marzari N, Novoselov KS, Geim AK, Ferrari AC: 
Uniaxial strain in graphene by Raman spectroscopy: # peak 
splitting, Griineisen parameters, and sample orientation. Phys Rev B 

2009, 79:205433. 

1 2. Ni ZH, Yu T, Lu YH, Wang YY, Feng YP, Shen ZX: Uniaxial strain on 
graphene: Raman spectroscopy study and band-gap opening. ACS 
Nano 2008, 2(1 1):2301 -2305. 

1 3. Mohr M, Papagelis K, Maultzsch J, Thomsen C: Two-dimensional 
electronic and vibrational band structure of uniaxially strained 
graphene from ab initio calculations. Phys RevB 2009, 80:20541 0. 

1 4. Lu Y, Guo J: Band gap of strained graphene nanoribbons. Nano Res 

2010, 3(3):189-199. 

1 5. Mei H, Yong Z, Hong-Bo Z: Effect of uniaxial strain on band gap of 
armchair-edge graphene nanoribbons. Chin Phys Lett 201 0, 
27(3):037302. 

1 6. Datta S: Quantum Transport : Atom to Transistor. Cambridge: Cambridge 
University Press; 2005. 

1 7. Ahmadi MT, Ismail R, Tan MLP, Arora VK: The ultimate ballistic drift 
velocity in carbon nanotubes.J Nanomaterials 2008, 
2008(2008):769250. 

1 8. Wong J-H, Wu B-R, Lin M-F: Strain effect on the electronic properties 
of single layer and bilayer graphene. J Phys Chem C 201 2, 
116(14):8271-8277. 

1 9. Liao WH, Zhou BH, Wang HY, Zhou GH: Electronic structures for 
armchair-edge graphene nanoribbons under a small uniaxial strain. 
Eur Phys J B 201 0, 76:463-467. 

20. Sun L, Li Q, Ren H, Su H, Shi QW, Yang J: Strain effect on electronic 
structures of graphene nanoribbons: A first-principles study. J Chem 
Phys 2008, 129(7):074704. 

21. Chang CP, Wu BR, Chen RB, Lin MF: Deformation effect on electronic 
and optical properties of nanographite ribbons. JAppI Phys 2007, 
101(6):063506. 

22. Huang M, Yan H, Heinz TF, Hone J: Probing strain-induced electronic 
structure change in graphene by raman spectroscopy. Nano Lett 
2010, 10(1 0):4074-4079. 

23. Shah R, Mohiuddin TMG, Singh RN: Giant reduction of charge carrier 
mobility in strained graphene. Mod Phys Lett 5 201 3, 27(03):1 350021 . 



doi:1 0.1 1 86/1 556-276X-8-479 

Cite this article as: Johari and Ismail: The effect of uniaxial strain on 
graphene nanoribbon carrier statistic. Nanoscale Research Letters 201 3 8:479- 



Submit your manuscript to a SpringerOpen 0 
journal and benefit from: 

► Convenient online submission 

► Rigorous peer review 

► Immediate publication on acceptance 

► Open access: articles freely available online 

► High visibility within the field 

► Retaining the copyright to your article 



Submit your next manuscript at ► springeropen.com 



v. 



J 



